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ABSTRACT 

We present a regularized maximum likelihood weak lensing reconstruction 
of the Deep Lens Survey F2 field (4 deg^). High signal-to- noise ratio peaks in 
our lensing significance map appear to be associated with possible projected fil- 
amentary structures. The largest apparent structure extends for over a degree 
in the field and has contributions from known optical clusters at three redshifts 
{z ~ 0.3, 0.43, 0.5). Noise in weak lensing reconstructions is known to potentially 
cause "false positives" ; we use Monte Carlo techniques to estimate the contami- 
nation in our sample, and find that 10-25% of the peaks are expected to be false 
detections. For significant lensing peaks we estimate the total signal-to-noise ra- 
tio of detection using a method that accounts for pixel-to-pixel correlations in our 
reconstruction. We also report the detection of a candidate relative underdensity 
in the F2 field with a total signal-to-noise ratio of ~ 5.5. 



Subject headings: gravitational lensing 
universe — galaxies: clusters: general 



dark matter — large-scale structure of 



1. Introduction 



The mass function o f galaxy clusters is a well known probe of cosmological parame- 
ters (IHaiman et al.ll200ll ). Since cluster rnass is not a direct observabl e other proxies suc h 
as optical galaxy richness flKoester et al.M2007l). X-ray te mperature (jRosati et al.l l2002l ). 
the Sunyaev-Zel'dovich decrement (jCarlstrom et. al.l 120021 ). or weak gravitational lensing 
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( ISchneiderlll996l ) are used to trace the mass distribution. To produce accurate constraints 
on cosmological parameters the mass-observable conversion must be properly calibrated in 
order to avoid potential biases. With weak lensing, mass overdensities are selected from maps 
of the projected mass density and peaks of mass are identified as a function of their signal- 
to-noise ratio. The observable in this case is the weak lensing 'shear' which is independent 
of baryonic physics, and theoretically well understood. This is a major strength that weak 
lensing has over other r nethods. Known limitations of this method a re projection effects due 
to large scale structure (IHennawi fc Spergelll2005l : IWhite et al.ll2002l ) which can produce spu- 
rious detections, a low peak detection significance, complicated point-spread function (PSF) 
corrections, and the observational expense (since it requires deep imaging). In spite of these 
limitations, forecasts for future imaging surveys predict that this method will still provide 
large samples of shear-selected galaxy clusters which will allow fo r constraints on cosmologi- 
cal parameters, including dark energy (IMarian &: Bernsteinll2006l ). In addition to individual 
clusters, studies have also pointed out that the global statistics of dark matter peaks in wea k 
lensing maps can be used to put constraints on cosmology (jjain fc Van Waerbekd I2OOOI ) . 
Within a given dataset cosmological information from peak statistics could also play a com- 
pleme ntary role in breaking de generacies in cosmological parameters from studies of cosmic 
shear (iGavazzi &: Soucaill 120071 ). 



The current generation of deep optical imaging su rveys are creatirig the fi rst ever maps 
of the dark matter distribution with weak lensing. In IWittman et al.l (120061 ) we presented 
preliminary maps from the Deep Lens Survey (DLS), covering a total of area of 8 deg^ spread 
over five fields. Initial maps from the Canada-Fr ance-Hawaii Telesc o pe Le gacy survey cover- 
ing four separate 1 deg^ fields were presented in iGavazzi fc Soucaill (120071 ) . Re cently, results 



from the Hubble Space Telescope COSMOS 2 deg^ field were reported in iMassey et al. 



(I2OO7I ). To reconstruct the projected mass density (k) in these surveys from weak lensing, 
all of these studies have implemente d a 'd irect' reconstruction method based on or related 
to the technique of iKaiser fc Squired (119931 ). This reconstruction method is computationally 
efficient but has a number of known limitations; for instance, it requires setting an arbitrary 
smoo thing scale, and it is difficult to incorporate other lensing effects such as magnifica- 
tion JSeitz et al.lll998l l Several different reconstruction techniques have been developed to 
overcome these issues, the most promising of which advocate reconstructing the deflection 
potential on a grid rather than directly determining k. These techniqu es have pr i marliy 
be en applied to pointe d observations of individual clusters, for instance in |Jee et al.l (120071 ) 



or 



Bradac et al.l (120061 ). In this study we use a regularized maximum likelihood technique to 



reconstruct the projected mass distribution over a wide area, four square degree field from 
the DLS. 



Our paper is organized as follows: In [|2] details of the imaging data, PSF correction. 
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and the source galaxies used in our analysis are presented. In ^ the maximum likelihood 
algorithm used in our reconstruction and its application to the DLS are discussed. In ^the 
structures observed in the lensing reconstruction are presented along with their likely optical 
counterparts. In ^ the expected false peak rate due to lensing shape noise is discussed as 
well as our method of determining the total signal-to- noise ratio of lensing peaks. In ^ our 
results are summarized and we discuss future directions of our work. 



2. Data 



The DLS is a deep imaging survey of five widely separated four square degree fields 
(F1-F5) in BVRz' . Fields in the DLS were selected in a blind manner the only restrictions 
were to avoid brig ht foreground galaxies , known low redshift clusters {z < 0.1), and areas 
of high extinction (jWittman et al.ll2002l ). A primary motivation for this was to create an 
unbiased shear selected sample of clusters. In this study we restrict our analysis to the F2 
field centered on R.A. = 09^19""32.4% decl. = -^30°00^00^^ Observations of the F2 field were 



carried out using the MOSAIC I imager (IMuUer et al.l 119981 ) on the Kitt Peak Mayall 4m 
telescope. Observations of F2 began at Kitt Peak in 1999 November and ended in 2004 
November. This was the first DLS field and to have complete imaging. The observing 
strategy for the DLS is to split each field into a 3 x 3 grid of 40' x 40' subfields. Each 
completed subfield consists of twenty 900s R band exposures, and twenty 600s exposures 
in each of the BVz' filters. The R band is used as the primary filter to measure galaxy 
shapes for weak lensing, and observations were only carried out in this filter when the seeing 
FWHM was < 0.9". The remaining filters {BVz') are primarily used to measure photometric 
redshifts of source galaxies, however these are not used in our current analysis. 

For our lensing analysis the shapes of background galaxies are measured from the co- 



added R band subfield 
given in 



rame s. A detailed discussion of the DLS weak lensing pipeline is 



Wittman et al.l (120061 ) here we provide a brief overview. Basic reductions such as flat 



fielding and bias are performed using the IRAF package MSCRED. The MOSAIC cameras 
each consist of eight individual CCD detectors, and the PSF in each detector is determined by 
identifying stars from the stellar locus in magnitude-size space. We fit the spatial variation of 
the PSF in each detector usi ng a third order polynom ial, and a 3 x 3 rounding kernel is used 
to circularize the PSF as in lFischer fc TysonI ( 19971 ). The PSF corre cted R band images o f 
each subfield are co-added using our custom software DLSCOMBINE ( IWittman et al.ll2006l ). 
The final co- added subfield images reach a depth of i? ~ 26. 



To initially detect objects in each subfield image we use the SExtractor (IBertin fc Arnouts 



19961 ) package. Shapes for detected objects are then remeasured with adaptive moments us- 
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ing ELLIPT O (ISmith et al.ll200ll : lHevmanal2006l ). which is a partial implementation of the 
algorithm of Bernstein fc JarvisI ( 2002 ). Objects which triggered error flags in ELLIPTO 
were removed since these are objects where adaptive moments did not converge or indicated 
a problem with shape measurement. We also rejected objects which triggered SExtractor 
flags indicating a fatal error occured during shape measurement, the object contained a satu- 
rated pixel, or the object was too close the edge of the image. For our analysis we use source 
galaxies in the magnitude range 22.0 < R < 25.5, where the number counts in the field peak 
at -R ~ 25.5. We include only galaxies with ELLIPTO size > 1.25 x PSF at the position of 
each galaxy. The maximum ELLIPTO size is set at 20.0 pixeP in order to eliminate poten- 
tial low redshift, low surface brightness galaxies from our sample. Our size and magnitude 
cuts produce a sample of source galaxies that each have a signal-to-noise ratio > 15. We 
also reje ct objects whose obse rved ellipticity is e > 0.6, as these are likely superpositions of 
objects ( iWittman et al.ll2000l ). The dilution (or seeing) correction is performed using simu- 
lations described in §3.2[ This is chosen over the analytic dilution correction described for 
the ELLITPO method in iHeymand (120061 ). since it has be en pointed out that th is correction 
does not correctly incorporate the kurtosis of the PSF (IHirata Sz SeljakI |2003| ). The final 
subfield catalogs were stitched together to create a supercatalog of source galaxies for the 
F2 field. This catalog contains ~ 328, 000 galaxies over the field, or ~ 23 galaxies arcmin^^. 



or 



Our current source galaxy catalo g is not the same catalog used in lWittman et al.l (120061 ) 



Khiabanian Sz DelF Antonio! (120081 ). In both of these studies the catalog was restricted to 



brighter galaxies {R < 25) with larger galaxy sizes than we use d here. This work uses a cat- 
alog w ith a factor of two more source galaxies than was used in iKhiabanian fc Dell' Antonio 



3. Weak Lensing Analysis 
3.1. Formalism 

— * 

In the thin lens approximation the mapping of light from the source plane P to image 
plane 9 is given by the lens equation (3 = 9 — a{9). The deflection angle a of light due to a 
two dimensional Newtonian deflection potential ip oi a, lens is given by 

a{9) = v^. (1) 

The surface mass density (k) of the lens can be calculated directly from the deflection 
potential by 

^ = 15(^^11+^22), (2) 

^crit ^ 
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where subscripts on if) refer to the partial derivatives ipij = 
mass density (Scrit) is given by 



-'crit 



Here the critical surface 



(3) 



where D^, D^, and Dds are the angular diameter distances to the lens, source, and between 
the lens and source respectively. The tidal gravitational field of the lens, or shear, is described 
in complex notion by 7 = 71 + 172 where the shear components are given by 



7i = Ui'u - V'22) 



(4) 



72 = ^12 = ip' 



'21 



and the complex reduced shear is given hj g = 'y/d — n) fjSchneider et al.ll200Cir ). For a more 



detailed discussion of lensing we refer the reader to iKochanek et al.l ( 20051 



3.2. Regularized Maximum-Likelihood Reconstruction 

Maximum l ikelihood methods for w eak lensing cluster mass reconstruction were first 



Bridle et al. 


(1998 


) and 


Seitz et al. 


(1998) 



To produce a weak lensing convergence map of the F 2 field we u s e a re gularized maximum 



likelihood ap proach based on the method develop ed by lSeitz et al.l (119981 ) with improvements 
described in iKhiabanian fc Dell'Antonid (120081 ). In this technique the deflection potential 



in the field is determined over a grid, and the convergence map is generated directly from 
the deflection potential map using equation ([2]). 

To produce a convergence map with dimensions N^xNy, a grid of the deflection potential 
with dimensions (2A^2; + 4) x [2Ny + 4:) is used. The twice larger grid is required here to fix the 
ringing effects in the projected mass maps caused by second order numerical differentiation 
of the deflection potential. The extra rows and columns are needed to compute k and 7 
using second order finite differencing. With this method we minimize the function J^, 



:f = -x^ + A7^, 



(5) 



where A is the regularization coefficient, and TZ is the regularization function. The term 
is determined from 

1 ^ (efc- < e > (xfc))^ 

aTZ. ^ ' (6) 



X 



k=l 
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where Ng is the number of source galaxies, ek is the complex ellipticity of a galaxy at 
position Xk, and < e > is the expected ellipticity at Xk- We compute the expected ellip- 
ticity distribution and the dispers ion 0"^, as a function of reduced shear from simulations 
(IKhiabanian &: DeirAntonidl2008l ). The ellipticity distribution i n our simulations is based 
on the measured distribution from the Hubble Ultra Deep Field (IBeckwith et al.ll2006l ). We 
regularize using a zeroth-order regularization function given by 



n 



m=l n=l 



Pmn) 



(7) 



where is the convergence at a grid point and pmn is the prior. This form of the regulariza- 
tion function is chosen for simplicity, but also ensures the smoothness of the reconstruction. 
The regularization coefficient. A, in J-" represents a compromise between the best fit and 
the closest m atch to the prior. We minimize the function T using the conjugate gradient 
method from lPress et al.l (119921 ). 



Our reconstruction proceeds at a series of different resolutions, beginning with a coarse 
20 X 20 grid of the potential and a completely uniform prior. This outputs a coarse potential 
which is used to create a convergence map. The resulting convergence map is smoothed 
and used as a prior to the next level of res olution. Using a smoothed pri or has been shown 
to produce more accurate r e constr uctions (ISeitz et al.lll998l : iLucy Ill994j ). As described in 
Khiabanian fc DelF Antonio! (120081 ) we choose the regularization coefficient at each of the 
higher resolution reconstructions by minimizing T (eq. E]) with multiple values of A (between 
and 10) along with minimizing only TZ (eq. [7]). We scale to values between and 1, 
using its lowest value obtained when A = and its highest value obtained when only TZ is 
minimized. Similarly, we scale TZ to values between and 1. The intersection of the scaled 
vs TZ curve and the scaled line = TZ determines the proper value of the regularization 
coefficient. This processes of producing maps at a given resolution and determining the 
proper regularization coefficient is repeated for three resolutions of the convergence map 
(20 X 20, 40 X 40, 80 x 80), until the final map resolution is achieved. Our final convergence 
map has dimensions of 80 x 80 with a plate scale of 1'.5 pixel"^. The resulting signal-to-noise 
ratio map (convergence map divided by the rms map) is shown in Figure [H The rms map 
(/^rms) is described further in §4.11 
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4. Dark Matter Peaks in F2 
4.1. K Signal- To-Noise Ratio Map 

We search for lensing peaks i n the F2 field by first c reating a k signal-to-noise ratio map 



(henceforth n — S/N map). As in iMiyazaki et al.l (120071 ) we construct a Krms map using 100 



Monte Carlo realizations of the original source galaxy catalog. We note here that half as 
many realizations could have been used and the following results would not have changed. 

In each realization, the position and ellipticity components of each galaxy are decou- 
pled and randomly assigned to new galaxies. Shuffling the source catalog in this man- 
ner is useful since it preserves the real variation in background galaxy source density. 
For each shuffled source catalog, a maximum likelihood k map is created using the same 
regularization coefficients at each resolution as was used to create our original k map 
khiabanian fc DelFAntoniol boosh . 7rom the resulting set of 100 Monte Carlo maps we 



created the Krms map shown in Figure [2l Each pixel in this map represents the la standard 
deviation at this point over the set of Monte Carlo reconstructions. An — S/N map for F2 
is made by dividing the original k map by the Krms map (Figure [T]) . Contours in the k — S/N 
map range from a peak signal-to-noise ratio (z/) of 1.0 — 6.0; negative contours have been 
omitted here for clarity. 



We detect peaks in the k — S/N map using the SExtractor package (IBertin fc Arnouts 



19961 ). Peaks are separated by setting the minimum contrast parameter in SExtractor to 
which separates all possible peaks. We measure peaks relative to a zero background level 
and use the height of each peak to approximate the signal-to-noise ratio of detection. We 
set the minimum area for detection at 9 contiguous pixels and the minimum threshold for 
detection at 0.01. Our definition of a peak is similar to that given in Jain & Van Waerbeke 
(200), however in our case our threshholding criteria eliminate many low-significance peaks 
based on size and significance. This use of a minimum threshold provides a cleaner peak 
catalog because peaks below this level of significance do not correspond to real detections. 

The edge pixels in our k — S/N map correspond to the edge of the F2 field, where the 
imaging is shallower and the noise is higher. Therefore, there are fewer galaxies that corre- 
spond to these pixels. Because the Monte Carlo procedure preserves the spatial distribution 
of galaxies, these pixels are systematically underconstrained. As a result the variance along 
the edge of /trms map is smaller than in the rest of the field. This can be seen as the dark 
border in Figure [21 We find this causes spurious detections along the edge of the k — S/N 
map, so to remove these we filter out a border region of size 3' from the lensing peak catalog. 



The resulting peak distribution in F2 as a function of u is shown in Figure [31 The 
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solid black histogram shows the distribution of positi ve peaks, the most signifi cant being the 
confirmed cluster Abell 781. Unlike previous studies (iGavazzi fc Soucailll2007l ) our use of an 
effective peak threshold cuts off the low signal-to-noise ratio end of the peak distribution, 
however below this level peaks do not correspond to real detections. The dashed histogram is 
the distribution of negative peaks where the k— S/N map has been flipped to the positive axis 
for comparison. The physi cal interpretation of s igniflcant negative peaks are underdensities 
in the matter distribution (IMiyazaki et al.ll2002l ). A possible underdensity is detected in the 
DLS F2 fleld, which we discuss further in §5.4[ 



4.2. High S/N Peaks 

To be consistent with pr evious studies we restric t our study of lensing peaks in the F2 



fleld to peaks with v > 3.5 (IGavazzi fc Soucaill 120071 ). Above this signal-to-noise ratio we 



detect 12 peaks, with positions and u values listed in Table [H All of these high signal-to- 
noise ratio peaks appear to lie along projected structures resembling fllaments in the k — S/N 
map shown in Figure [H Positions of peaks are overlaid in circles on the n — S/N map, with 
the rank of each peak (sorted by peak value) shown in each circle. Three distinct projected 
structures appear in the reconstruction — the largest projected structure spanning over a 
degree. Below we discuss each of these structures and their association with likely optical 
counterparts in the imaging data. 



4-2.1. Eastern Structure 



In the Eastern section of our map we detect a projected structure consisting of eight 
signiflcant peaks which extend along the North-South direction and span for over a degree. 
The two highest signiflcance peaks (Peaks 1 & 2) are associate d with the cluster com plex 
Abell 781, of which we reporte d a previous lensin g detection in IWittman et al.l (120061 ). X- 
ray observations of this system (jSehgal et al.ll2008l ) have shown this complex consists of four 
separate clusters. In our lensing reconstruction we have separated two distinct clumps, where 
the easternmost clump in our map (Peak 1) is the most signiflcant {u = 6.6). Due to the 
pixel scale in our reconstruction (1'.5 pixel"^) this clump is detected as a superposition of 
two clusters, CXOU J092110+302751 at z = 0.42 7, and CXOU J09205 3+302800 at ^ = 
0.291. Our previous direct lensing reconstruction in IWittman et al. J2006h referred to these 
as clumps C and B respectively. To the West of this peak we detect our second highest 
ranked peak (Peak 2 with u = 5.7) which is t he cluster CXOU J09 2 026+302938 a t z = 
0.302. This peak was referred to as clump A in IWittman et al. J2006h . ISehgal et al.l J2008h 
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recently reported the X-ray detection of an additional cluster to the west of the A781 complex 
(XMMU J091035+303155 at z = 0.428) which has a lensing significan ce of la - 2a when 
fit to a Navarro, Frenk, & White (NFW) profile (INavarro et al.lll996l ). In our maximum 
likelihood reconstruction we note that we are not able to detect this cluster in our map. 
As we discuss in §5.31 more work is needed to explain th e low detection s i gnific ance of this 
cluster. Images of the Abell 781 complex can be found in lWittman et al.l (120061 ). 



South of the Abell 781 complex along the Eastern structure we detect a new significant 
peak (Peak 3) with v = 5.3. The DLS optical imaging reveals a large excess of galaxies 
located near the peak center. Photometric reds hifts of likely clust er members taken from 
the SDSS Data Release 6 (DR6) photoz2 table (lOyaizu et al.l 120081 ) place this cluster at a 
redshift of z ~ 0.5. We stress that this is only a preliminary redshift based on a small 
number (~ 5) of likely cluster members. This lensing peak likely has contributions from 
other cl usters along the line of sight; for instance, we note that the SDSS Maxbcg cluster 
catalog (IKoester et al.l 120071 ) which overlaps the F2 field also detects another cluster with 
optical galaxy richness Ng^h = 15 and redshift z = 0.297 located ~ 3'. 8 to the southeast of 
this peak position. At this angular separation however, this cluster does not likely contribute 
much signal to this lensing peak. 

South of this peak we detect two other significant peaks in the Eastern structure. Peaks 
7 and 12, with u = 4.1 and u = 3.6 respectively. The DLS optical imaging reveals several 
galaxy overdensities in the vicinity of peak which likely contribute to the lensing signal in 
each. We currently do not have richness or redshift information for the groups, so cannot 
disentangle the relative contributions to each lensing peak. South of this grouping we detect 
another grouping of three peaks : Peaks 4 (z/ = 4.5), Peak 6 {u = 4.2), and Peak 9 {u = 3.9), 
where the fourth ranked peak is north of the two other lensing peaks. From the DLS imaging 
Peaks 4 and 6 do not appear to correspond to galaxy overdensities and are potentially false 
positives. Peak 9 is likely associated with many optical overdensities; for instance, ~ 1.5' 
to the southwest of peak 9 there is a group scale systen i which has a single spectroscop ic 
measurement from the SDSS DR6 of z = 0.339 ± 0.0002 jAdelman-McCarthv et aP 120081 ). 



The Eastern Structure appears to be a superposition of structures lying at several differ- 
ent redshifts [z ~ 0.3,0.43,0.5). More follow-up work is needed however to disentangle the 
relative strengths of the le nsing contributions at each redshift. The Smithsonian Hectospec 
Lensing Survey (SHELS) (iGeller et. al.ll2005l ). a redshift survey being conducted in the F2 
field, should provide further insight into this structure. 
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4-2.2. South- West Structure 



In the southwest portion of our map we detect another structure containing three sig- 
nificant peaks. Two of these peaks (Peaks 8 fc 11) are associated w ith the system DLSCL 
J0920. 1+3029 which we previously reported in IWittman et aLl (120061 ). This system consists 
of two known clusters separated by ~ 9', which each have a spectroscopic redshift of z = 0.53. 
The northern peak (11th ranked peak) is detected with u = 3.7 and the southern peak (8th 
ranked peak) is de tected with u = 3.9. I mages of the clusters associated with peaks 8 and 
11 can be found in IWittman et al.l (120061 ). 



North of these peaks we detect a new peak (Peak 10) with u = 3.8. This lensing 
peak lies near an overdensity of galaxies ~ 3' to the East centered on R.A. = 09^^16^15. 7*^, 
decl = +29°50'05.5". A ~ 3' separation between the center of the lensing peak and the 
optical overdensity is large and could possibly indicate that this peak is spurious, however 
the broad lensing contours do encompass this galaxy overde nsity. Photon i etric redshifts of 
likely cluster members from the SDSS DR6 photoz2 table (lOyaizu et al.l |2008| ) place this 
cluster at a redshift of 2; ~ 0.55. This redshift is also only an estimate based on a small 
number (~ 5) of likely cluster member galaxies. This cluster and the other two lensing peaks 
likely place the redshift of much the southwest structure at z ^ 0.5, however other structures 
along the line of sight could also be contributing to the lensing signal here. 



4.2.3. Central Peak 



Near the center of our map we detect a new significant lensing peak (Peak 5) with 
u = 4.4. The optical imaging reveals two brightest cluster galaxies (BCG's) coincident with 
this peak, surrounded by an excess of fainter red galaxies (Figure H]). The BCG located at 
R.A. = 09^18'^36M, de cl = +29°53^08^^ has a spectros copic redshift from the SDSS DR6 
of z = 0.3171 ± 0.0002 jAdelman-McCarthv et al1l2008l ). We note that the SDSS Maxbcg 
cluster catalog detects a cluster ~ 2'. 5 to the East of this peak with richness A^gais = 17 at a 
slightly lower redshift z = 0.278, which could also contribute to the lensing signal. 



4.3. Low S/N Peaks 

We comment here on lower S/N peaks which appear to be associated with the projected 
structure seen in the lensing reconstruction, but do not fall in our high S/N sample. South 
of our 5th ranked peak, the k — S/N map (Figured]) reveals a structure which we detect as 
two peaks; the first located at R.A. = 09^19"'26.1^ decl. = +29°22'06.4" with u = 2.7 and 
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the second located at R.A. = 09''19'"19.0" decl. = +29°33'04.2" with v = 2.4. 

North of Peak 5 we detect another structure in the n — S/N map which are resolved as 
three separate peaks, which we designate as A, B, & C. Peak A is at R.A. = 09^20™56.5'^, 
decl. = +30°16'36.6" with u = 2.7, Peak B at R.A. = 09^20"^55.0^ decl. = +30°08'44.0" 
with u = 2.7, and Peak C at R.A. = 09^20"'02.9^ decl. = +30°06'34.6" with u = 3.4. The 
superpositions of these clusters as well as potentially other clusters along the line of sight 
likely contribute to the lensing peaks in this structure. 



5. Discussion 



5.1. False Peak Rate 



The noise in weak gravitational lensing reconstructions is influenced by several sources: 
the distribution in background galaxy positions, ellipticities, and redshifts, spatial variations 
in the sky background and PSF, and pixel-to-pixel correlations due to the reconstruction 
methods. As a result, the noise distribution in the lensing maps is strictly non-Gaussian, and 
there is the potential for "false positives" — peaks in the lensing map that do not correspond 
to any real objects. 

Several previous studi es have addres s ed the issue of the false peak rate in weak lens- 
ing mass reconstructions. iHamana et al.l ( 120041 ) used N-body simulations and ray-tracing 
to construct catalogs of expected clusters to compare with them. They reported that the 
completeness and contamination of the catalogs depended on the S/N threshold, and that 
even at relatively high threshol ds (~ 5a), a 20 - 30% incompleteness and a similar contami- 
nation ra te would be exp ected. ISchirmer et al.l (120071 1 used a variation of the aperture mass 
statistic (ISchneiderlll996l ) and estimates of the projected galaxy overdensity to estimate the 
noise contamination rate in a sample of ESO-WFI fields. When they varied the S/N thresh- 
old, th ey found contamina t ion ra tes of 30-80%, depending on the depth and seeing of the 
fields. iHetterscheidt et al.l (120051 ) applied two variations of the aperture mass statistic to 



50 deep Very Large Telescope fields; they set their thresholds so as to approximate a 20% 
contamination rate from noise peaks in their sample. Even though these studies used very 
different reconstruction methods, data sets and selection criteria, they all concur in finding 
20-40% contamination in the rough S/N range that we consider. At lower S/N, the very large 
number of "detections" in both k and randomize d map s indicate a much higher contami- 
nation rate. We also mention that I Van Waerbekd (120001 ) studied the noise in weak lensing 
reconstructions, however their analysis is not directly applicable to us since it does not apply 
to maximum likelihood reconstructions which use a regularization term. 
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T o estimate the po tential false peak rate in our high S/N sample due to lensing shape 



noise (iKubo et al.ll2007l ) we performed two tests. In our first test we created a mean map of 
the 100 Monte Carlo reconstructions described in §4.11 Each pixel in this map is the average 
value of this pixel over the individual 100 Monte Carlo map realizations. An average signal- 
to- noise ratio map is created by dividing this map by the mean of the Krms map (/trms/lO) 
from §4.1[ Peaks are detected in this average signal-to-noise ratio map using exactly the 
same procedure used to analyze the original k — S/N map. The resulting peak histogram is 
shown in Figure E] (Left). We find that above our high signal-to-noise ratio cut {u > 3.5) 
three peaks are recovered. In a second test we created individual signal-to-noise ratio maps 
from each of the 100 Monte Carlo realizations. Here each individual Monte Carlo map is 
divided by the K^ms map from §4.1[ Peaks are detected in each of these individual signal-to- 
noise ratio maps using the same detection procedure described in §4.11 The resulting average 
histogram is shown in Figure (Right). Here we detect ~ 1 false peaks above our high S/N 
cut. Both of the above tests indicate a potential false peak rate due to shape noise of 1-3 (or 
~ 10 — 25%) above our high signal-to-noise ratio cut {u > 3.5). This percentage compares 
reasonably well with an estimate based on the proximity of our peaks to overdensities of 
galaxies in physical and redshift space. Based on the analysis of the environment near the 
peaks, we see 2-3 (16-25%) of the structures appear not associated with or at > 3' from 
identified galaxy concentrations. 

In the above tests we do not simulate the signal due to uncorrected large scale structure 
along the line of sight which will increase the number of detected objects that are falsely 
interpreted as cluster candidate. The result of these tests therefore are only an estimate to 
the lower bound of the false peak rate. We emphasize that although false peaks are expected 
above our high signal-to-noise ratio cut, the majority of our peaks appear to be associated 
with real clusters which also appear to trace the large-scale filamentary structure. 

We n ote here that res i duals in the PSF correction could also potentially cause false 



peaks. In iKaiser fc Squired (119931 ) type convergence maps this can be probed at some level 



by creating a B -mode map where so urce galaxies are rotated by 45 degrees relative to each 



grid point (e.g. iMassey et al.l (120071 )). One limitation of our maximum likelihood algorithm 
however is that in its current implementation it cannot create a B-mode map. To probe PSF 
systematics we have alternatively measured radial shear B-modes for our high significance 
peaks (Figure [6]). In the figure we show profiles for four of these peaks spread over the 
structures in the F2 field. In each panel the profile is center on the peak coordinates given 
in Table [TJ The measured B-mode for each peak is consistent with zero, indicating that 
there is no residual contamination coming our PSF correction. Profiles for the other high 
significance peaks in the field are similarly consistent with a zero B-mode signal. 
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5.2. Total S/N 

To be consistent with previous studies we have selected peaks from the k — S/N map as 
described in §4.11 This provides a good method to globally select peaks from the map, how- 
ever because our reconstruction technique does not track pixel-to-pixel correlations, neither 
the peak nor the integrated signal-to-noise ratio are good estimates of total signal-to-noise 
ratio of detection. To calculate the total signal-to-noise ratio for peaks we use a different 
approach where the original k map and the individual Monte Carlo maps are block averaged 
at different levels of block size, k — S/N maps at each block size are then created in the same 
manner used to create the original k — S/N map ( §4.11) . For small values of block size, due 
to the pixel-to-pixel correlations in the n and Krms maps, both the signal and noise increase 
linearly, hence the peak signal-to- noise ratio (z/) will stay relatively constant. When the 
block size becomes larger than the noise correlation length, noise will increase proportional 
to square root of the block size; therefore, the peak signal-to-noise ratio will grow as more k 
pixels are used in the average. At some level of block size the peak signal-to-noise ratio will 
be maximized when the pixel size matches the characteristic size of the cluster. The peak 
signal-to-noise ratio at this block size can then be used to estimate the total signal-to-noise 
ratio of detection. Past the maximum as the block size continues to increase we expect 
that the signal-to-noise ratio will decrease as more noise pixels are added. An example of 
this technique is shown in Figure [71 here we show peak S/N vs. block size curves for three 
different peaks in our map. The top curve is for the peak grouping 1 & 2, the middle curve 
is for peak 5, and the bottom curve is for the peak grouping of 8, 10, & 11. In each of 
the curves the peak S/N reaches a maximum at certain level of block size, and the value 
of the peak S/N at that block size is used to approximate the total signal-to-noise ratio of 
detection. This provides a near-optimal measure of the integrated signal-to-noise ratio for 
isolated clusters, subject to the constraint of finite pixel sizes. One limitation of our total 
signal-to-noise ratio estimation is that for large block sizes some peaks can become blended 
with other peaks. This behavior is exhibited in the top curve (Peaks 1 & 2) and the bottom 
curve (Peaks 8, 10, & 11). In such cases we can only estimate the total signal-to-noise ratio 
for the combined system. We list total S/N values using this method for the remaining peaks 
in Table [H 

Our method of estimating the total signal-to-noise ratio for lensing peaks is qualitatively 
similar to previously-used techniques in this field. For example, the P-statistics introduced 
by Schirmer et al. (2007) use variations in the scale of the kernel to effectively identify peaks 
based on their signal at a smoothing scale matching the peak size, even though they do not 
take the additional step of identifying the S/N at thi s kernel with the total S/N of the peak . 



Our method is also similar to the method used in iKaiser. Squires, fc BroadhurstI (119951 ): 



however in that paper, the method was used to estimate the S/N ratio of faint galaxies 
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rather than in estimating the significance of the lensing signal itself. 



5.3. The Abell 781 Complex 



As mentioned in §4.21 there is a known X-ray cluster associated with the western portion 
of Abell 781 complex that does not appear in our reconstruction. Our lensing map does hint 
at a peak at this position but this peak appears to occupy only 1-2 pixels in our map. We 
recent ly reported a Icr — 2cr lensing detection of this cluster from an NFW fit in lSehgal et al. 



( 120081 ). The lensing signific ance of this cluster is low compared to what is expected from the 
X-ray mass of this cluster (jSehgal et al.l 120081 ). The low lensing significance could possibly 
be due to a remaining systematic, or could potentially point toward an interesting physical 
effect which causes X-ray and lensing masses to be discrepant. Observations of this region 
with a different telescope/imager are planned, and a detailed analysis will be presented in a 
future paper. 



5.4. A Candidate Void? 



In §4.1l we mentioned the detection of a candidate underdensity in the k— S/N map with a 
significant negative signal-to-noise ratio. Strictly speaking, this region is a relative underden- 
sity since it is underdense relative to other structure in the map. Our abih ty to detect a rela- 
tive u ndersity is unaffected by the mass sheet degeneracy transformation (jSchneider fc Seitz 
1995h . The undersity is located at R.A. = 09^20'^22.7^ decl. = +29°54'03.5" and is detected 
with a peak significance ~ 5 and a total signal-to-noise ratio of ~ 5.5. From our lensing 
reconstruction alone we cannot directly confirm that this underdensity is in fact a void. A 
galaxy redshift survey, an alternative ma pping of structure, is a well known method of di- 
rectly tracing underdense regions or voids ( Geller &: Huchra 1989[). W e are conducting such a 
survey in the F2 field, the SHELS redshift survey (IGeller et. al.ll2005l ). which should provide 
additional insights into this underdense region. Confirmation of an underdensity detected in 
a lensing reconstruction would demonstrate that lensing maps can trace mass underdensi- 
ties, in addition to oy erdensities. We note that an excess of negative peaks was reported in 
Miyazaki et al.l (120021 ). however this excess has not been confirmed with a redshift survey. 
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6. Summciry 

In this work we have presented the apphcation of a regularized maximum hkehhood weak 
lensing reconstruction to the Deep Lens Survey F2 field. Our lensing reconstruction reveals 
several projected structures in the F2 field which appear to lie along a network of potential 
filaments. In the southwest section of our map we detect a structure which appears to lie 
at a redshift z ~ 0.5. Another structure, the Eastern-most structure in our reconstruction, 
contains the confirmed cluster Abell 781 and appears to extend southward for over a degree. 
Known optical clusters at three redshifts z ~ 0.3, 0.43, 0.5 are potentially associated with 
significant lensing peaks along this structure. Near the center of our map we detect another 
significant peak at a redshift z ~ 0.32. The cluster associated with this lensing peak lies 
at a roughly the same redshift as many of the peaks in the Eastern structure, however it 
is unclear if this structure is physically associated with any part of the Eastern structure. 
These two structures also appear to surround a relative underdense region with i/ ~ 5, and 
total signal-to- noise ratio ~ 5.5. 

The noise in our reconstruction is estimated from 100 individual Monte Carlo realiza- 
tions of the source galaxy catalog. We find that generating Monte Carlo k maps with this 
reconstruction technique is computationally expensive, but does allow us to create a statis- 
tically correct k — S/N map. Monte Carlo realizations also aid in estimating the number 
of expected false peaks due to lensing shape noise. We performed two tests here using the 
Monte Carlo k maps; in one test we created an average Monte Carlo significance map, and in 
another test we created individual Monte Carlo significance maps and averaged the results. 
Both of these tests estimate 10 — 25% false positive detections due to shape noise in our high 
peak signal-to-noise ratio sample. This is comparable to the number of peaks in our dataset 
which lack close counterparts in galaxy overdensities. 

We have also presented a new method to estimate the total signal-to-noise ratio of 
detection for lensing peaks. Because pixel-to-pixel correlations are not traced in our lensing 
reconstruction (or in any traditional direct reconstruction) the peak signal-to-noise ratio 
underestimates the total lensing significance. Our method of estimating the total lensing 
significance uses block averaging of the k map and the individual Monte Carlo maps to 
create k — S/N maps with different levels of block size. The peak signal-to-noise ratio will 
become maximized here when the peak matches the characteristic size of the cluster, and 
this peak signal-to-noise ratio can be used to estimate the total lensing significance. We find 
that this method works quite well for isolated peaks, but for a grouping of peaks only the 
significance of the combined group can be measured. 

In future work, photometric redshifts of source galaxies in the F2 field could allow us 
to make tomographic maps using this method and would also allow for the determination 



- 16 - 



of the mass of significant lensing peaks. Magnification information in tlie F2 field could 
also be incorporated into our reconstruction, potentially breaking the mass sheet degeneracy 



( iBroadhurst et al.lll995l ). We are also pursing optical cluster finding in the DLS, however 



this is beyond the scope of our current study. 
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Table 1. Lensing peaks with z/ > 3.5 in the DLS F2 field 



Peak Rank 


K.A. ( ) 


Uecl. [ j 




z 


b/I\t 


1 


140.2294 


+30.4677 


6.6 


0.427, 0.291 


7.0 


2 


140.0829 


+30.5070 


5.7 


0.302 


7.0 


3 


140.3072 


+30.2299 


5.3 


~ 0.5 


5.4 


4 


140.4161 


+29.6915 


4.5 


N/A= 


4.0 


5 


139.6471 


+29.8726 


4.4 


0.317 


4.9 


6 


140.2994 


+29.5624 


4.2 


N/A= 


4.0 


7 


140.5350 


+30.1093 


4.1 




3.6 


8 


138.9965 


+29.5569 


3.9 


0.53 


3.9 


9 


140.4174 


+29.4782 


3.9 


~ 0.34 


4.0 


10 


139.0007 


+29.8431 


3.8 


0.55 


3.9 


11 


138.9588 


+29.6458 


3.7 


0.53 


3.9 


12 


140.5499 


+29.9367 


3.6 




3.6 


void 


140.0946 


+29.9010 


- 5 




~ 5.5 



'^u is the peak signal-to- noise ratio measured from the n — S/N 
map 

^S/Nt is the total signal-to- noise ratio measured using the 
method outlined in §5.21 

'^Peaks 4 and 6 do not appear to be associated with galaxy 
overdensities and are likely false positives 

Note. — Due to the resolution of our map, some peaks are 
blended together. In this case we report the total signal-to-noise 
ratio of detection for the combined systems. 
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Fig. 1.— Weak lensing k - S/N map of the F2 field centered on R.A. = 09^19™32.4^ 
decl. = +30°00'00". Postive contours of signal-to- noise ratio (white) are overlaid and enclose 
the 4, 2, 3, and la regions; negative contours (red) enclose the —4, —3, —2, and —la regions. 
Black circles indicate the positions of the most significant peaks (z/ > 3.5) with the rank of 
each peak shown at the center of each circle. 
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Fig. 2. — The Krms noise map based on 100 Monte Carlo realizations of the F2 field. Each 
individual Monte Carlo map is created by shuffling the source galaxy catalog positions and 
ellipticity components. A maximum likelihood k map is made for each individual Monte 
Carlo realization. Each pixel in the K^ms map represents the la error at this point over the 
set of Monte Carlo reconstructions. 



-22- 



15 



10 - 












4 
V 



8 



Fig. 3. — Distribution of peak signal-to-noise ratio {v) measured from our S/N map of the 
DLS F2 field. The solid histogram is the distribution of positive signal-to-noise ratio peaks, 
the dashed histogram is the distribution of negative signal-to-noise ratio peaks. Negative 
values have been flipped here to the positive axis. 
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Fig. 4. — Galaxy cluster coincident with Peak 5 in our high signal-to-noise ratio peak sample. 
North is up, East is to the left. The cluster candidate contains two bright elliptical galaxies, 
the eastern-most having a spectroscopic redshift of z = 0.3171 from the SDSS DR6. 
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Fig. 5. — Distribution of peak signal-to- noise ratio {v) for the mean Monte Carlo map (Left) 
and for the average of the individual Monte Carlo maps (Right). As in our real lensing map 
solid histograms are positive peaks, dashed histograms are negative peaks (where negative 
values have been flipped to the positive axis). For our high signal-to- noise ratio cut (y > 3.5) 
three peaks are recovered in the mean map, ~ 1 is recovered in the average of the individual 
Monte Carlo maps. These tests provide an estimate of the number expected false peaks in 
our high signal-to-noise ratio sample due to lensing shape noise. 
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Fig. 6. — Radial profile of the B-mode component centered on four high-signficance peaks in 
our convergence map. The B-modes are all consistent with zero signal, indicating that the 
residuals in the PSF correction are not introducing a spurious signal near each lensing peak. 
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Fig. 7. — Examples of total signal-to-noise ratio estimation using the method outlined in 
§51 The top curve is for the peak grouping 1 & 2, the middle curve is for peak 5, and the 
bottom curve is for the peak grouping 8, 10, & 11. In each, the peak signal-to-noise ratio 
{u) is measured for different levels of block size after block averaging the k and individual 
Monte Carlo maps, k — S/N maps are then created for each level of block size. The curve 
for Peak 5, an isolated peak in our map, exhibits a clear maximum. For the two other 
curves the maximum is not as well defined due to the proximity of nearby peaks, however an 
estimate for each peak grouping can still be obtained. As expected the total signal-to-noise 
ratio values are higher than our estimates using only the peak signal-to-noise ratio from our 
K — S/N map. 



